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Abstract 

We use boundary field theory to describe the phases accessible to a tetrahedral qubit coupled to 
Josephson junction chains acting as Tomonaga-Luttinger liquid leads. We prove that, in a pertinent 
range of the fabrication and control parameters, an attractive finite coupling fixed point emerges 
due to the geometry of the composite Josephson junction network. We show that this new stable 
phase is characterized by the emergence of a quantum doublet which is robust not only against 
the noise in the extemal control parameters (magnetic flux, gate voltage) but also against the 
decoherence induced by the coupling of the tetrahedral qubit with the superconducting leads. We 
provide protocols allowing to read and to manipulate the state of the emerging quantum doublet 
and argue that a tetrahedral Josephson junction network operating near the new finite coupling 
fixed point may be fabricated with today' s technologies. 
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1. Introduction 

Quantum impurity models [1] provide a natural paradigm to describe a large number of non- 
perturbative phenomena occurring in one dimensional quantum devices such as point contacts, 
constrictions, crossed quantum wires and Josephson junction (JJ) chains [2, 3, 4, 5, 6]. While 
a standard perturbative approach is accurate when the impurity is weakly coupled to the environ- 
mental modes needed to fully describe the quantum device, there are situations in which impurities 
are strongly coupled to such environmental modes: when this happens, it is impossible to disen- 
tangle the impurity from the rest of the system, the perturbative approach breaks down and one 
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has to resort to the non perturbative tools provided by boundary field theories (BFT) [1,7], which 
have been shown to yield accurate descriptions of many realistic low dimensional systems [8]. 

Prototypical non perturbative impurities states are realized in systems exhibiting the Kondo 
effect [9, 10, 11, 12], or in situations where static defects appear in Tomonaga-Luttinger liquids 
(TLL)s; for both settings, a renormalization group approach leads, after bosonization [8], to the 
emergence of boundary sine-Gordon models[13, 14]. In both cases, the interaction with the impu- 
rity makes the boundary coupling strength to scale to a stable strongly coupled fixed point (SFP) 
which is characterized by a fully screened spin in Kondo systems or by the eff"ective disappearance 
of the impurity in a "healed" TLL [15, 16]. More remarkable states are achieved when a finite cou- 
pling fixed point (FFP) - characterized by new non-trivial universal indices - emerges; for instance, 
this happens in the overscreened Kondo problems [9, 10, 11, 12], or in crossed TLLs where, as 
a result of the crossing, some operators turn from irrelevant to marginal, leading to correlation 
functions exhibiting power-law decays with nonuniversal exponents [4, 5, 17]. 

Quantum impurities are realizable also in superconducting Josephson devices [18, 19, 20, 21]. 
Superconducting Josephson chains with a weak link [18, 22] and SQUIDs [19, 23] may be indeed 
described by boundary sine Gordon models yielding a phase diagram with only two fixed points: 
an unstable weakly coupled fixed point (WFP), and a stable one at strong coupling. At variance, 
for pertinent architectures of the Josephson junction network (JJN) one may find a range of fabri- 
cation and control parameters where a stable FFP emerges in the phase diagram [6] allowing for 
the engineering of superconducting devices exhibiting enhanced quantum coherence [24] and 4e 
superconductivity [20]. The stable fixed point is associated with the emergence of a doubly de- 
generate ground state which may be regarded as a quantum doublet described by a spin 1 /2 degree 
of freedom coupled to the plasmon modes via the boundary interaction; as a result, one may use 
superconducting devices not only as good candidates for the design of solid state quantum bits 
[25] but also as efficient quantum simulators of the various physical behaviors realizable in Kondo 
systems. 

Engineering quantum doublets robust against noise and decoherence is of paramount impor- 
tance for applications to quantum information processing. Quite recently a remarkably robust 
two-level quantum system has been shown to emerge in a device made with six JJs arranged (see 
Fig.l) in a symmetric tetrahedral geometry [26, 27]. In Refs.[26, 27] it has been pointed out that, 
when each intemal loop is pierced by a dimensionless magnetic flux f = n, the ground state is 
doubly degenerate and that the degeneracy is robust against small variations in the applied gate 
voltages and/or in the applied magnetic flux, since the degeneracy is split only to second order in 
the charge and flux noise. Remarkably, the design of a tetrahedral quantum bit [26, 27] may be 
modified so as to make it robust also against the noise in the Josephson junction energy [28]. 

For any practical purpose (control or reading out the state of the quantum bit) one needs to 
connect the tetrahedral quantum bit to external leads, which may induce new decoherence effects 
spoiling the coherence of the quantum doublet [29]. As we shall show, realizing the leads with 
TLLs, enhances the coherence of this device. This happens if one uses one-dimensional super- 
conducting leads which may be mapped onto TLLs with the Luttinger parameter g depending on 
the fabrication parameters of the JJN [18, 19, 22, 23, 30, 31]. In this paper, we address this issue 
by analysing the device made by the central tetrahedral JJN - T - depicted in Fig.l connected to 
three JJ chains ending in bulk superconductors at fixed phase j = 1, 2, 3; the resulting network 
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is depicted in Fig. 2. We show explicitly that the interaction of the central region T with the low- 
energy collective excitations (plasmons) of the leads merely renormalizes the "bare" parameters 
of T and, thus, does not break the tetrahedral "symmetry" responsible for the robust groundstate 
degeneracy [26, 27]. As a result we exhibit a JJN which is robust also against fluctuations of the 
lead's parameters and argue that such network may be fabricated with nowadays technolgies [32]. 
The paper is organized as follows. 

In section 2, we derive an efi'ective spin-1/2 Hamiltonian for T (Fig.l). In agreement with the 
results of refs.[26, 27], we find that the spectrum of the efi'ective Hamiltonian admits a twofold 
degenerate groundstate and argue that, in the effective theory, T may be regarded as a spin-1/2 
degree of freedom So- Finally, we derive the low-energy description of the JJN depicted in Fig. 2. 
We show then that the efi'ective theory is a l-i-l dimensional field theory, where the central region 
is described by a spin interacting with the leads via a pertinent boundary interaction. 

In section 3 we use a perturbative approach to account for the couplings between T and the 
TLL leads realized with JJ chains. Here, we use a renormalization group approach to determine 
the flow of the running boundary coupling strengths, as a function of the system size and we 
determine the range of values of g and / where the perturbative approach breaks down. 

In section 4, we analyze the large-scale behavior of the JJN depicted in Fig. 2 in the strong 
coupling limit where the boundary efi'ects induced by T become relevant. In particular, we find 
the set of minima of the boundary interaction potential at the SFP and construct the instanton 
operators connecting two degenerate minima. Using the renormalization group approach [33], we 
analyse the flow of the coupling constants associated to the instantons. Finally, we show that - for 
~ 10 ^ / ~ ^ ~ fo ^ < g < 3 - a finite coupling fixed point emerges. 

In section 5 we show that instantons are responsible for the emergence of a two-level quantum 
system which, due to the network's architecture, is robust, for -y^ ^ / ~ ^ ^ fo' against the 
decoherence arising from coupling the cental region to the leads. Moreover, we show how the 
state of the emerging quantum doublet may be manipulated by acting upon the extemal control 
parameters. 

In section 6 we point out that the quantum doublet emerging in a tetrahedral JJN operating 
near the FFP may be realized with today's technologies. 
Section 7 is devoted to concluding remarks. 

The appendices provide the necessary mathematical background to follow the analysis carried 
in the text. 

2. Boundary field theory description of the tetrahedral JJN 

In this section, we derive the low-energy, long-wavelength description of the JJN depicted 
in Fig. 2. We shall find that this JJN may be described by a l-l-l dimensional field theory, with 
a pertinent boundary interaction describing the central region. In the following we shall firstly 
derive the Hamiltonian associated to the central region T and then construct the BFT describing 
the JJN in which T is connected to TLL leads realized with JJ chains ending with three bulk 
superconductors held at fixed phases (pi,(p2, (P3- 
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Figure 1 : The central region T. The junctions are assumed to have the same Josephson energy Ey^ while the three 
loops are threaded by the dimensionless fluxes f\ , f2, /a, respectively. 




Figure 2: The tetrahedral JJN. The central region T, depicted in Fig.l, is made with four superconducting grains 
connected by six quantum Josephson junctions tuned nearby the degeneracy between two charge eigenstates. T is 
connected, via a coupling A «: Ew, to three leads realized with one-dimensional Josephson junction arrays fabricated 
with junctions of nominal Josephson energy Ej and ending in three bulk superconductors held at fixed phases (p\,ip2, 



2.1. The central region T 

T is fabricated with six quantum Josephson junction, joined to each other as depicted in Fig.l. 
We assume that, for each junction, the charging energy Ec is much bigger than the Josephson 
energy E^/^. To prevent Coulomb blockade from forbidding charge transport across T, we further 
assume that to each superconducting grain is applied a gate voltage Vg, tuned nearby the degener- 
acy between the charge eigenstates with charge equal to N and N + I. Under these assumptions 
one can describe each superconducting grain with a quantum spin-1/2 variable Sq^ [18, 19, 21]; 
this leads to an effective spin-1/2 representation for the Hamiltonian "Kt describing the central 
region. 

If each internal loop in Fig.l is pierced by a (dimensionless) magnetic flux / (i.e. /i = /i = 
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f-i = /), a standard procedure [21] yields 



3 ,^ r 3 



i=0 ^ 



i=l 



(1) 



(=1 

In Eq.(l), Sq\ i = 1, 2, 3, denotes an effective spin-1/2 operator lying on the three outer sites of T, 
while Sf^ lies on the central site. The parameter H oc e*Vg - N - ^ is determined by the uniform 
gate voltage bias. 

"Kt can be exactly diagonalized: its eigenvalues and eigenstates are given in Appendix A. 
There we show also that, for | < / < y and \H\ Ew, the ground state is twofold degenerate; 
here, we denote the states of the degenerate quantum doublet as |0, 1) = Ifl") and |0, 2) = The 
effective spin-1/2 operator Sg, acting onto the two-dimensional subspace spanned by {If]'), |JJ,)} may 
be then represented as 

% = \Z\^)(^'\r''y ^ (2) 

(T,(T 

where a = 1 ... 3, r are the Pauli matrices and cr ='[\, JJ.. 

Although a spin-1/2 degree of freedom emerges quite naturally in pertinently engineered JJN 
[6, 20, 21], its robustness - against any detuning of e*Vg off the degeneracy value A'^ -I- ^ as well 
as against any small deviation of the dimensionless magnetic flux / from its optimal value n - 
is a very challenging task. Indeed in all the devices analyzed in Refs.[6, 20, 21] the twofold 
degeneracy is realized only if / is fine-tuned to n, since any displacement from f = n breaks the 
doublet degeneracy already to the first order in /-;r. At variance, due to the frustration induced by 
the presence of the central spin SqO), the tetrahedral central region T is much more stable against 
noise in the external control parameters. 

2.2. The emerging quantum doublet 

A simple symmetry argument provides us with an hint on why the architecture of the central 
region is relevant for the robustness of the emerging quantum doublet. For this porpouse we first 
prove that this robustness depends crucially on if the number of junctions needed to fabricate the 
central region is even or odd. 

If the central region is realized with an odd number of junctions -as it happens, for instance, 
in the Y-shaped networks analyzed in [6]- the ground state degeneracy is a consequence of the 
invariance of the effective Hamiltonian under time-reversal symmetry T which, for a three-spin 
effective Hamiltonian describing a triangle-shaped central region [6], is explicitly realized as 

3 

f = {Y[^j\Tn, , (3) 



5 



where the operator [Y\%i Cy] (^"^ being the first Pauli matrix) changes the sign of each one of the 
three spins, while fi23 reverses the label order of the spins (that is, (123) — > (321) '). Since, 
when the number of spins is odd, the ground state has total spin 1/2, and since T reverses the 
sign of the total spin, one sees immediately that the ground state must be twofold degenerate. 
Such a degeneracy is, however, easily spoiled by a displacement of the applied gate voltage off 
the degeneracy point, i.e., by a nonzero value of the parameter H, which is coupled linearly to the 
total spin of the state: any detuning of e*Vg off the degenerate value A'^ + ^ breaks then the ground 
state degeneracy already to the first order in e*Vg - N - ^. At variance, when the central region 
is effectively described by a spin Hamiltonian with an even number of spins (i.e., the number of 
junctions needed to fabricate the central region is even) the ground state is a spin singlet and, thus, 
insensitive to a nonzero value of H. Although robustness versus accidental displacements in the 
applied gate voltage is guaranteed, to get a degenerate ground state requires a fine tuning of the 
applied flux f to f = n: any displacement of the applied flux off the optimal value f - n breaks 
the degeneracy between the two ground states already to the first order in / - ;r. 

The tetrahedral geometry of the central region "Kx provides an optimal compromise between 
the above complementary issues, since the twofold ground state degeneracy is protected by sym- 
metries that are realized for any value of /. Indeed, time-reversal T is now realized as = IsT^s 
where I5 reverses the sign of all the spins and T123 is defined in Eq.(3). Looking at the explicit 
form of the states | tt) J ID^ given in appendix Appendix A, one may readily see that 

Ill) = l5l^) , I^) = l5lll) , (4) 

and that 

fu3\^\) = e-'^\^\) , ri23iu) = e^iii) . (5) 

Since I5 is not a symmetry of 'Hx, one concludes that I5I H) must be an eigenstate of 'Hx, inde- 
pendent of I 11), but degenerate in energy with it, as explicitly shown in Appendix A. This ensures 
that the twofold ground state degeneracy is allowed within a rather large range of values for the 
applied flux /. 

Robustness against variations in the external flux may be easily seen also with the help of a 
perturbative analysis. Indeed, in order to break the degeneracy between | ff), | JJ.) one may look at 
situation where the three fluxes are split according to fj ^ f + Sj (with \Sj/f\ «c 1). Writing 
'Hj[{f + 6i}] as 

E ^ 

■HtH/ + 6t}] « 7Yx[{/}] - i^e'f J][Si{S^^ns];^\] + h.c. + 0[{6^}] , (6) 

^ ,-=1 

one may easily compute matrix elements of '7Yx[{/ + Si}] - 'Hriif}] within the subspace spanned 



'Notice that multiplication of [Yl]=i cr]] by 7'i23 is needed in order for the symmetry to be preserved also when a 
magnetic flux is applied 
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by I It), I JJ.), provided one knows that: 



{cr\ J'^^ ^ sin(/) J] 5, (7a) 
1=1 1=1 

3 

{a\ S^'V]I<t) oc sin(/) (e'f (5i -62 + e^'^^) , (7b) 



r=l 



with cr ='[]■, U, and 11 =U,, |], =11. From Eqs.(7), one sees that, for f = n, the ground state stays 
degenerate at least up to first order in the displacements {6j} 6j. To second order in the {6j}, one 
finds that the low energy eigenstates | t|') and | JJ.') and their corresponding eigenvalues are given 
by: 

e^> = -EM + ^-^ ^ -A , (8) 



where 



_ 1 2 i 



^{6] - 61) + {Si - 6]) + (61 - 6]) 



A = ff V2 ^{61 - 6\? + {6\ - 6lf + (61 - 6^)^ . (9) 

From Eq.(8,9) one sees that the second order corrections to 'H[{f}] may be recasted into an effec- 
tive Hamiltonian 'H±, describing the central region as a transverse magnetic field term, given by 



•H^ = 5^cos(^)S^ + 5^sin(^)S-^ , (B^ = -A) . (10) 

Of course, when all the 6j are equal then 5^ = 0- 

In the next section we shall derive the BFT describing the central region T connected to three 
TLL leads realized with three JJ chains fabricated with junctions of nominal Josephson energy Ej 
and ending in three bulk superconductors held at fixed phases (pi,(p2, (fs- 

2.3. Connecting T to TLL leads: the boundary Hamiltonian. 

To derive the BFT describing the tetrahedral JJN depicted in Fig. 2 we require at first that the 
leads are realized by one-dimensional JJN for which the Josephson energy Ej is bigger than the 
charging energy Ec; we further assume that there is an uniform gate voltage Vg acting on each 
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junction tuned nearby the degeneracy point between the states with N, or N + I Cooper pairs at 
each junction. As a result [18], each lead may be described by a one-dimentional spin- 1/2 chain; 
if each chain is made out of L sites, the Hamiltonian describing the leads is given by 



Leads 



L-1 



£1=1,2,3 y=o,L-i 



a= 1,2,3 j=0 



with 



. d 



CT' 



a,] 



-I 



(11) 

(12) 
(13) 



Here, (pa.j is the is the phase of the superconducting order parameter at site- j of the a-chain, Pg is 
the projector onto the subspace of the Hilbert space with either N,or N + I Cooper pairs at each 
superconducting grain, and is the effective strength of the charge interaction between nearest- 
neighboring junctions. As a result the low-energy long wavelength limit of the Hamiltonian (11) 
can be described by a one-dimensional spinless TLL Hamiltonian, given by 



= ±y r 



dx 



1 (d(^a\ (d^a 
+ U 



(14) 



u\dt) \ dx 

where describes the collective plasmon modes of the leads, while the Luttinger parameters g 

(A = (E- - ?>E]l\6Ec)IEj, a is the lattice 



and u are given by ^ = 



— A—, u — uE I 

2(ff-arccos(^))' •' 



L\2 



2 arccos(^) 

step) [18, 22]. 

The leads are connected to T by means of three Josephson junctions, of nominal strength 
A «c Ew < Ej, connecting the endpoints of the leads to the outer sites of T. The Hamiltonian 
describing this interaction is given by [6] 



3 

-K, = -Aj]{(Sfre-—+iS^^re'—) 

i=l 



(15) 



Using Eq.(15) together with the spectrum of T given in Appendix A, the Schrieffer- Wolff (SW) 
procedure yields an effective boundary Hamiltonian involving only the low-energy degrees of 
freedom of T. A rather lengthy computation yields 



74 =2£i/2cos 

+ 4£3S^X 
+ 4E,%J] 



O,(0)-O,^i(0) 
V2 



-I- 



cos 



0,(0) - O,-^i(0) n 
V2 ^2, 



-I- 



cos 



sm 



'2 


2) 




3^(J- 


cos 


2 ,. 


2) 




3^(7- 


COS 



OX0)-Oy+l(0) 

V2 

Oy(0)-Oy+l(0)' 

V2 



-I- 



BnSr 



(16) 
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where 

E.= ,f'' . ,E,= ^/f' . ,E, = 2E, , (17) 

and 

_ 24V3^W-^) .... 
Bii ^2 • ^^^■> 

One sees from (16, 17, 18) that the term fiycr^ explicitly breaks - to the second order in the control 
parameters g and /- the degeneracy between | H) and | JJ,). In the following sections we shall show 
how the low-energy plasmon modes of the leads renormalize the parameters of "Hb. 

3. Perturbative analysis near the WFP 

In this section, we determine the flow of the running boundary coupling strengths, and argue 
about the emergence of nonperturbative fixed point in the phase diagram accessible to the tetrahe- 
dral JJN. 

3.1. Renormalization group flow of the boundary coupling near the WFP. 

To check the stability of WFP, we derive the renormalization group (RG) equations for the 
running boundary coupling strengths. To do so, we use of a boundary version of the RG approach 
to perturbed conformal field theories, developed by Cardy [34]. The starting point is given by the 
Euclidean boundary action corresponding to the Hamiltonain in Eq.(16), which is given by 

dr cos (dj ■ i^(r)) + 4E, dr % '^"^ («i " + 2 ) + 
4^3 dr '^"^ (^^0' - 2) j cos ((?y • ^(r)) + 4E, dr J] (^^0' - 2) j cos (dj ■ ;^(t)) , 



+ 



with /3 = ikBTr\ 



(19) 



^^=•0) ' ^^ = (vi/2) ' ^^ = (-75/2) • ^^^^ 



We have defined a^(t) = xiO,T) = [Xiix,T),X2ix,T)], with^i(;c,T) = ^[Oi(;c,r) - 02(;c,t)], 
X2(x,t) = -^[Oi(x, r) + 02(jc, r) - 203(x, r)], in order to evidence that the "center of mass" 
field 0(jc, r) = ■^[^iix,T) + 02(x, r) + 03(x, r)] decouples from Sg^ as expected from charge 
conservation at T [4, 5, 6]. One may then write the "free" action Sq = ^LeadD?] + Ss[^, O] only in 
term of the fields and;tf2- Namely, 



g 

'LeadDr] = ^ J ^"^ J 



uXdr 1 \dx 



(21) 
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and 



JO(t), 



1 



5s[O,0] = -- J Jr-^[l-cos(0(T))]-- J JT{5||COs(0(T))+5^sin(0(T))cos(O(T)-^)} , 

(22) 

being the imaginary time action for the quantum spin variable So reported in appendix Appendix C 
(0, are the polar angles: see appendix Appendix C for details). The partition function for the 
JJN, Z., is then given by 



(0) 



where 



^ ;=1,2 ^ 



SLeadl;?]-5s[<I>,0] 



(23) 



(24) 



denotes the imaginary time time-ordered product, and the boundary interaction action is written 

as 



,^(T) . (25) 

(?eA a=x,y,z 

In Eq.(25) the colons : . . . : denote normal ordering with respect to the vacuum of the bosonic 
theory [18], with all the expectation values computed with respect to the WFP Hamiltonian, and 
A = {+«!, +a2! it^s}. Furthermore, Eq.(25) defines the dimensionless couplings g"^. 
Expanding Z in the couplings yields 



Z 



^ (feA a,b=x,y,7. 

-J[ Z YJ^^^/i III drdr'df 



;(r))+ 



(26) 



ajj,yeA 



" { : : : e'^'^^^^'^ : : e''^'^^^'^ : )<S;^(t")S^(t')S^(t)) + 



(27) 



where we have used the fact that ( : : . . . : ; ) is different from only if Xi"=i (^^^^ = 0. 

Following the procedure of Refs.[19, 33], one finds that a rescaling of L implies a renormalization 
of the couplings, according to the RG flow equations that, to second order in the boundary coupling 
strenghts, are determined by the short-distance operator product expansions (O.RE.s) 



u\t-t'\ 

z 



(28) 



when di + dj + dk = 0, and 



u{t - t') 



u{t - t') 



u dr 



(29) 



10 



as well as by the O.P.E.s between the spin- 1/2 operators which, may be derived starting from the 
equations of motion for a spin in magnetic field B = + Bl + B^, given by: 

'Bt\ 



S^(T) = S^(0) 
+ S-^(O) 

- S^(0) 
S-^(T) = S-^(O) 

+ %iO) 
-S^(0) 

+ S^(0) 

- S^(0) 



2 cos^ (p sin^ 9 sinh^ {'y) '^°^^(^''')) 
i cos 6 sinh(5T) - sin(20) sin^ 6 sinh^ (~2~)) 
/ sin sin sinh(5T) + sin 20 cos sinh^ ( ~2~ )) ' 
-2 sin^ (f> sin^ 0sinh^ {~^) ^c<^sh(5T)j + 

/ cos sin 6 sinh(5T) - sin (f> sin(20) sinh^ (~2~)) 
/ cos 6 sinh(5T) + sin(2^) sin^ 6 sinh^ (^)) ' 



1 +2sin2 sinh^ (—'!) + 



|/ sin sin 9 sinh(5T) - cos (p sin(20) sinh^ ( ~2~ )) 
^« cos 4> sin 9 sinh(5T) + sin (p sin(26') sinh^ (^)) ' 
where the angles and 9 are defined by 

B\\=B cos 9, B^ = B sm9cos(p , By = Bsm9s'm(p . 
The resulting O.P.E.s are: 

{S^(t')S^(t)),,^. 



|S-^(t')S-^(t) 
|S^(t')S^(t) 



S^(t')S^(t) 
S^(T')S^(r) 



S^(r')S-^(T) 



1 B 

- + -(r - t) [S^(t) cos 9 + SI(t) sin sin 9] 
1 B 

- + -(r' - r) [S^(t) cos + %{t) cos sin 



4 2 
1 5 

4^2 

^S^(t) + ^(t'-t) 



J. 

- + — (r' - t) [S^(t) cos sin + S^(t) sin (/> sin 
B 



2 

2' 
B 

2' 



2Sg(t) + t(t'-t) 
^Sg(t) + ^(t'-t) 



-I cos 9 + S^(t) sin sin 
-I sin cos (p + S(j(t) cos 9 



-I sin sin + S^(t) cos sin 



-^S-^(r) + |(T'-T) 



--I sin sin + S^(t) cos 9 
i 

--I cos 9 + S;^(t) cos sin 9 



B 

1' 



jSi(r) + -(r' - T) 



--I COS sin + S^(t) sin sin 9 



(30) 



(31) 



(32) 
(33) 



(34a) 
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When 5|| = Bj_ = spin operators do not evolve in imaginary time and Eq.(34a) reduces to the 
usual spin algebra 

S^S^ = i^"^! + i/6"^^S^ . (34b) 

Since at the WFP the fields x{t) and Sg(t) are decoupled, from Eq.(28,29) one finds that two 
diff"erent contributions to the boundary action are generated to the second order in the boundary 
coupling strenghts; to evidence these contributions, one rewrites Eq.(19) as 



J Jo 

with the So-dependent couplings SjiSc] and £^[80] given by 



(35) 



GjiSc] = Ei+ liE^G + 2^3 |cos 
fijLSc] = El - 2iE^S}j + 2E3 \ cos 



2n 

yO-2) 

2n 

yO-2) 



+ sin 



Sq + sin 



2n 

y(i-2) 

2n 

-jij-2) 



(36) 



£;[Sg] and £^[80] contribute to the boundary action in Eq.(25) through the O.P.E. as 



^fiy[SG(r)] : e-'"r^^^^ : fi,v[SG(r')] : e-''r>l^^'^ 



-AiE,Ei%{T) - 2E1E3 cos 



2n 



S'q(t) + sin 



: e''^"-'^^'^ :} , (37) 



|(*-2, 



S^T) 



and 



J]{Bj[Sg(t)] : e'''r^^^ : 6j[Sg(t')] : e-'^r^^^'^ : +6j[Sg(t)] : e-'^r^^'^ : fi,[SG(T')] : e'-^^'^^^^') :} 



m(t - r') 



cos 



|u-2) 



- sin 



|0-2, 



ZQf; 



M dr 



(38) 



From Eq.(38), one sees that, though terms proportional to |^ are not present in Eq.(19), they are 
dynamically generated by the renormalization group procedure, yielding a new term given by 



6S 



(I) 



dr <cos 



2n 

-0-2) 



SUt) - sin 



2n 

-0--2) 



3 J 



u dr 



■ (39) 



The boundary coupling A is renormalized from the terms in Eq.(38). 
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For a set of (dimensionless) boundary couplings g^, the RG equations to the second order in 
the couplings are given by: 



dl 



(I -K) 81(1) + J] J] 

j8,7 b,c=x,y,z 



(40) 



where / = ln(L/Lo), Lq is a reference length scale, is the scaling dimension of the boundary 
operator associated to the coupling and the coefficient C^'^'^ are defined by the O.P.E.s. Thus, 
deriving the coefficients C^'!^'^ from Eqs. (37,38) and taking into account that h = ^(h = l)for any 
vertex operator (for ^^p), one finds that the RG equations for the dimensionless 

couplings Gi(L) = L^~~^Eu G,(L) = L^~'^E-, G^{L) = L^~~^E^, Ga(L) = A, are given by 



JGi(Z) 

dl 
dG-Sl) 

dl 
dG,(l) 

dl 
dG^m 

dl 



1 

1 - - 



1 

1 - - 



Gi(/) - 2 



Gl(l) 



MGuG„G,,Ga)) 



GS) + 4Gi(/)G,(/) + 2G3(/)Ga(/) = MGuG,, G,, Ga) 
G,(l) + 2Gi(/)G3(/) + 2G,(/)Ga(Z) = /3,(GuG„ G,, Ga) 
= 4a(/)G3(/) = MGuG,, Gs, Ga) • 



1 

1 - - 



(41) 



From Eqs. (41) one sees that no linear term appears in the jS-function for the running coupling Ga, 
since the corresponding boundary operator has scaling dimension 1, while the ones containing 
the vertex operators have scaling dimension 1 /g. Furthemore, looking at the linear term in the 
jS-functions in Eqs. (41), one sees that the couplings Gi(/), G,(/), G3(/) are all irrelevant for g < 1 



and that Ga(0 scales roughly like e*^^"*^', so that, for g < 1, the term oc is irrelevant, as well. 
At variance, for g > 1, Gi(/), G^{l),Gj(l) are all relevant, while Ga(0 is only marginally relevant. 
As a result, is subleading with respect to Sg\ and, for g > I, one may write Eqs. (41) by 
keeping only leading contributions to the boundary interaction; thus one gets 



dGijl) 

dl 
dG,(l) 

dl 

dGAD 
dl 



i-i 



i-i 



1 

1 - - 



Gi(/)-2 



G\{1) 

G\{1) - Gl{l) - 



= /3i(GuG„G,) 



G,{1) + 4Gi(/)G,(/) = P,{GuG„ G3) 
G3(/) + 2Gi(/)G3(/) = yS3(Gi, G„ G3) 



(42) 



From Eqs. (42), one sees that, when g > \, the boundary interaction provides a relevant perturba- 
tion and that the pertubative approach breaks down when the leads reach the healing length , for 

i-i 

which the biggest dimensionless coupling is of order Lq Ej. Lq is the reference length such that 
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1-i 

Gi(l = 0) = Lq Ei, with the £,'8 defined in Eq.(17); since Eq.(17) shows that the biggest coupling 
is Gt,{1), Lh. should be determined from the requirement that G3 [ln(|^jj ~ ^ Ej yielding 

U - (^V ' Lo . (43) 



It should be noticed that, for g = 1, one should use Eqs.(41), instead of Eqs.(42), as all the 
four couplings Gi(/), GJJ), G^il), GfJJ) correspond to marginally relevant perturbations. To spell 
out the RG flows for g = 1, one should notice that G^(Lq) = if // = 0; from Eqs.(41), then, 
the RG flows run along the manifold in parameter space defined by G^{1) = Ga(/) = and by 
Gsil) = 2G\{1). Since G3 depends on / as 

— ^ = VG,{l)f , (44) 
dl 

one sees that the RG equations in Eq.(41) coincide with the ones obtained for the boundary cou- 
pling in the isotropic Kondo model. Thus, for g = 1 and G-(Lo) = 0, the tetrahedral JJN simulates 

a Kondo spin with isotropic couplings to the electrons and coincides with the Kondo length Lk, 

1 

given by ~ Loe^'s"^"'. A small value of H should just slightly change the RG flow, possibly to 
the one corresponding to a Kondo spin with anisotropic couplings to the band electrons [35]. 

In the next section, starting from Eqs.(41) we characterize the SEP and investigate its stability 
against the leading boundary perturbation. 



4. The Strong Coupled Fixed Point 

Erom Eq.(41) one sees that all the boundary couplings become relevant when g > \ \ this im- 
plies (43) that, as soon as L > , the running boundary couplings cross over towards strong 
coupling. At the SEP, the plasmon fields ;^(0) satisfy Dirichlet boundary conditions since their 
values must coincide with a minimum of "Kb; in addition, the leading boundary interaction is a 
combination of phase slip (instanton) operators describing tunnelling events between the neigh- 
boring minima of "TYb- 

4.1. Minima ofH^ 

Dirichlet boundary conditions are set by requiring that ;^i(0),;^^2(0) take values corresponding 
to a minimum of 'Hr. To determine the set of minima of "Kb it is most convenient to represent 
[4, 5] the spin- 1/2 operator Sg as a 2x2 matrix and rewrite 'Hr as 

^ ^(Vr(xm + VMO)) + Bn{f) VMO))-iVmO)) \ .... 

In Eq.(45) 5||(/) is given by Eq.(18) while, by comparison with Eq.(16), the V((y(0))'s are deter- 
mined as 
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j 

j 

j 

y,[;^ = -2£,^sin[(?,.-;^ . 



cos[a^- ■)(] 
cos[a^- ■)(] 



(46) 



One may now easily look for the eigenvalues of 'Kb as a function of ;^i(0) and ^2(0) obtaining 

AiCr(0)) = Vi - ^y? + y2 + (K + 5||)2 (47) 

A2Ct(0)) = Vi + ^y2 + y2 + (K + 5||)2 . (48) 
Since, for any value of one has that Ai < A2 the minima of "Kb are obtained from the minima 










*rwl 

a) 



Figure 3: Lattice of the minima of TYb. In the box: the triangle whose vertices may be made degenerate with a 
pertinent choice of the external phases, and the instanton trajectories associated to the operators . 



of Ai(;^(0)) only. We find that, in the;t'i(0),^2(0)-plane, these minima are located at the vertices 
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of the Kagome lattice displayed in Fig. 3. The polarization of Sg at each minimum is determined 
by minimizing the boundary energy after the V/s have been evaluated on the pertinent vertex of 
the Kagome lattice of minima. 

Since at the outer boundary (x = L) x\ ,X2 are connected to three bulk superconductors at fixed 
phases, the corresponding boundary conditions are given by ;k^i(L, t) = iXi,X2{l-, t) = fJ-2, with the 
relative phases /ii,/i2 being defined as /ii = (tpi - (fj)/ V2, and //2 = (^i + ^2 - '^'(pj)/ V6. Upon 
imposing Dirichlet boundary conditions for r) at both boundaries the mode expansion of the 
plasmon fields is given by 

;,,(x,r)=f,.^{(i-.)£p,-|]sin(?H)^^,-i..j , (49) 

where = XjiQ ^^id PjS are the zero mode operators. 

The eigenvalues of Pj are fixed by the boundary conditions since ^^Pj = they are 

given by 

(PuP2)i = ^f2g^-^ + ni2 + €i,-^ + ^{2ni3 - ni2 + 6i)j , (50) 

ni2, nu being relative integers. The index I accounts for the "red" (R), the "green" (G), and the 
"blue" (B) sublattices in Fig.3; €r = I, 6g = 1/2, 6b = 1/2; 6r = 0,6g = -1/2, 6b = 1/2, while the 
constants ^1,^2 cancel the terms oc //i,yU2 in;^i(0),;^2(0), respectively. The spin polarization of Sg 
is uniform throughout each one of the three sublattices; thus, to the minima of 'Hr are associated 
the states 



l^) = l 


\cos- 


|G) = | 


\cos- 


15) = 1 


i 

COS- 



4 



(51a) 
(51b) 
(51c) 



where cos 4r = . " and sin = ^ To the first order in (i.e. to the 



yi6+B2-B| 

^16+b2-B|I yi6+e2 2 ^16+Bj|-B|| yi6+B2 

first order in (/ - tt)) one has: 

"/-^^ (52) 

From the zero-mode spectrum, one sees that the energy of each field configuration gets a zero- 
mode contribution which is quadratic in p\ inserting the solution of Eq.(49) in the noninteracting 
Hamiltonian 'Hll, one finds the zero-mode contribution to the total energy for each of the three 
sublattices. Namely, one gets: 



• Red-points: 
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Zero-mode eigenvalues: 



(Pi,P2)r= ^/2gl^-^+nu + l,-^ + ^(2nl^-nu)j , (53) 



Zero-mode contribution to the total energy: 



17(0) _ 



nug 



-f^+n, + l 
2n 



+ 



-:r- + — p(2ni3 - nn) 
In V3 



(54) 



• Green-points: 

Zero-mode eigenvalues: 



13 -"12 - -j) 



(55) 



Zero-mode contribution to the total energy: 

nug 



^G,Mi2,ni3(/f) 



^1 i\' I H2 1 i;' 



(56) 



• Blue-points: 

Zero-mode eigenvalues: 



Zero-mode contribution to the total energy: 



EB,nu,ni-iiP) 



nug 



^1 ^ ^ ^2 ^ 1 ,^ ^ 1 



(58) 



Through the external phases /ii,)U2 one adds an effective quadratic potential breaking the degen- 
eracy of the zero-modes. It is quite remarkable that the zero-mode energy spectrum reported in 
Eqs. (54,56,58) displays, for convenient choices of the applied phases fJ.i,fi2, a threefold degener- 
acy between the three sites lying at the vertices of a single triangle (box in Fig.3). As we shall 
see in later sections, this feature of a tetrahedral JJN is crucial for its applications to quantum 
information processing tasks. 
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4.2. Leading boundary perturbation at the SFP 

In order to describe effects of the leading boundary perturbation at the SFP one should rewrite 
the boundary Hamiltonian in terms of the phase-slip operators representing instanton trajectories 
connecting nearest neighboring sites on the lattice in Fig. 3. For this purpose, it is most conve- 
nient to introduce the phase slip field operators @\{x, t), 02(^, t), dual to;^^i(;c, t),X2{x, t), and, thus, 
defined by [6] 



1 d&i{x,t) dxi{x,t) 



gdxi{x,t) d@i{x,t) 
u 



(59) 



ug dt dx u dt dx 

In particular, Eqs.(59) imply that &i 2{x,t) obey Neumann boundary conditions at jc = and at 
;c = L. As a result, the mode expansion of the instanton fields is given by 



®j{x, t)= ^|Tg\e^^ + ^Pj + i cos 



[nnx 



n*Q 



I L 



n I 



(60) 



with 



[PhOl^] = -iSij , [a'„,aiJ = n6n+m,oSij . (61) 

A quantum tunnelling between two adjacent minima lying on the lattice shown in Fig.3 involves a 
quantum jump between different zero mode eigenstates. The quantum phase slip operators corre- 
sponding to the allowed quantum jumps, Wi, Wj (i =1,2, 3), are given by: 



wl = 


:exp 




^(0) 


wl = 


:exp 


_ V3 


^(0) 


wl = 


:exp 




0(0) 



(62) 



The vectors Pi are defined so that -yf^g-^ is the "distance" between nearest neighboring eigenval- 
ues of P, as determined by the commutation relations 



[P, Wi] = 




=^fftWi , [P,W]] = 




which yield 



Pi 



V3/2 \ 

1/2 ; 



P2 





-1 



P3 



-V3/2 \ 

1/2 ; 



(63) 



(64) 



In Fig.3 we have represented the quantum jumps between the eigenvales of P corresponding to the 
operators Wj; the action of the hermitean conjugate operators Wi may be simply represented by 
reversing the arrows. 
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Since a tunnelling event between minima of "TYb involves a rotation in the two dimensional 
space spanned by the eigenstates of So, and since the allowed directions of tunnelling from a 
given minimum depend on the position of this minimum in the Kagome lattice, each quantum 
phase slip operator has to be multiplied by the spin operator mapping the state (i = R, G, B,) 
onto the other two states. As a result the leading boundary perturbation at the SFP is given by 

= -Ye-^Sl^Wl + Sl^Wl+Sl^Wl] + h.c. , (65) 

with Sij the operator sending the "j" spin state into the "i" one, while the parameters Y,y are 
computed in Appendix D. Since the scaling dimension of the operators {Wi, Wj} is | the running 
coupling strength for the dual boundary coupling may be defined as ^(L) = y{L)e'^ , with y{L) = 
LY(L). 

Even if ^ is, in general, a complex coupling strength, the renormalization group equations may 
be derived following the standard procedure used in [6] which, starting from Eq.(65), allows to 
determine the euclidean dual boundary action at the SFP as 

5b = -Ye'^ dr {SI^{t)WI + sI^(t)WI + 5;^(t)WJ} + h.c. , (66) 
Jo 

Here Wjir) = e with 0(t) = 0(0, /r); one may then compute the yS-functions for the 

boundary coupling strengths [6, 34] by resorting to the O.P.E.s 

u(t - t ) 



[WlS^]iT)[WlS^]iT') 



mS](.T) , (67) 



plus cyclic permutations. 

From Eq.(67), the second-order renormalization group equation for the running coupling ^(L) 
may be written as 

dm 



din 4^ 



^ = (l-|)OT-2e-^'W) , 



(68) 



where Lq is, again, a reference length scale. Eq.(68) is equivalent to the system of real differential 
equations for the real parameters y{L), y(L) given by 

dyiL) 



din 



U-=(l-^)y(L)-2cos(3y)y\L) 



yiD-^^^ = 2 sinOy{L))y\L) . (69) 

Mi) 

From Eqs.(69) one sees that the phase y is renormalized only to the second order in the boundary 
couplings and that, for y = kn/3 (k integer), there are lines of fixed points in the y - y-plane; 
furthermore, the line y = | is made of attractive fixed points. 

The phase diagram accessible to the tetrahedral JJN may be inferred from the RG equations 
near the WFP, derived in section 3, and the ones near the SFP, Eqs.(69). As displayed in Fig.4, 
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Figure 4: Phase diagram in the strip in the /-plane corresponding to < / < Itt: the weakly-coupled phase 
(corresponding to Neumann-(N)- boundary conditions at the inner boundary) is stable for g < I and for any /; the 
strongly-coupled phase (corresponding to Dirichlet-(D)- boundary conditions at the inner boundary) is stable for g > I 
and for any /. For g < 3 and - ■fj)<(/-^)^fga novel phase opens, corresponding to a stable finite coupling fixed 
point (FFP). The phase diagram for / i [0, In] is obtained by periodical extension of this picture. 

for 1 < ^ < 3 and -y^ ^ (/ - ^) ^ Y5 neither the WFP, or the SFP, is infrared stable. As a 
consequence, in this window of values of / and g, the infrared behavior of the JJN will be driven 
by an emerging FFP, whose properties will be analyzed in the next section. 

5. The quantum doublet at the FFP 

In this section we show that the renormalization of the instanton tunnelling strength Y, due to 
the interaction with the plasmon modes of the TLL leads, enforces the emergence of a quantum 
doublet robust not only against the noise in the external control parameters but also against the 
decoherence induced by the coupling with the plasmon modes in the leads. 

5.1. The emerging doublet at the FFP 

Since a real device has a finite size (L) the infinite degeneracy induced at the SFP by the 
eigenvalues of the zero-mode operators, is removed by finite-size effects; i.e., by the zero-mode 
contributions to the total energy. As a function of the external phases fl, the zero-mode energy 
associated to each eigenvalue p is given by 

<G,.K«.,..3[/^] = ^[^>i2,ni3)]^ . (70) 

From Eqs. (54,56,58), one sees that, for a pertinent choice of the phases ft, the three zero-mode 
energies associated to the vertices of a triangle as the one drawn in the box of Fig. 3, may be 
made degenerate. To explicitly show this, let us set /i* = y,jU2 = 0- For such a choice of the 



^This interval is determined in Sec. 6 and accounts for the fact that, for f n, the renormalization of Bj. may lift 
the degeneracy between the minima. 
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external phases, one gets Ef^Qlft] = E'^^qUT] = E'-^^^lft]. Thus, if one restricts himself to the 
three-dimensional subspace T of the Hilbert space spanned by the states \R), \G), \B), an effective 
low-energy description of the tetrahedral JJN at the SFP may be provided by the 3x3 Hamiltonian 
matrix "K^^ given by 



-Ye'^ -Ye-'^ 

-ye-'> -Ye'y 

-Ye'^ -Ye-'^ 



(71) 



Form (71) one sees that, when Y 0, the degeneracy between \R), \G) and \B) is broken 
and that the spectrum of the low-lying energy states admits as its groundstate a quantum doublet 
confined away from a singlet state. To see this, let us consider the two operators acting on T 
whose matrix representations are given by 



(72) 



While Ai corresponds to a cyclic shift of the states \R), |G), \B), A2 realizes the mirror inversion 
about the triangle height passing through the vertex corresponding to the |7?)-state. Though, as 
expected, when y = 7t/3, both Ai and A2 commute with "Kr, ['Hr,Ai] = {'Hr^Aj] = 0, they do 
not commute with each other since [Ai,A2] 0. This is enough to ensure that the spectrum of 
must contain at least one degenerate eigenvalue. To explicitly check it out, let us consider the 
eigenvalue equation associated to 'Hj^. It reads 








1 


1 




1 








Ai = 


1 





1 


, A2 = 












1 


1 
















-E^ + 3Y^E - 2 cos(3r)y^ = . 



(73) 



For 7 = n/3 there is a twofold degenerate eigenvalue E = -Y, and a non degenerate eigenvalue 
E = 2Y. The corresponding eigenstates are 

I - y)i = -^[\R) + \G) + m 

V3 

|2F) = \\R) + e'^G) + e~'^B)] . (74) 
V3 ^ J 

To rewrite using the states given in Eq.(74), one has to rotate it to Hr = P^'HrP, with the 
matrix P given by 



1 



1 1 



1 



1 e~'3" e'^" 
1 e''^" 



(75) 
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One obtains 



9L 9 



Similarly, by rotating A2, one gets 



-Y 

-y 








27 



(76) 



1 

1 
1 



(77) 



thus, A2 is the operator swapping the degenerate states \-Y)i,\- 7)2 with each other. 

Eq.(76) shows that a quantum doublet- this time robust also against decoherence induced by 
its coupling with the plasmon fields in the TLL leads- emerges in a tetrahedral JJN operated near 

the FFP. Indeed, accounting for the fluctuations of the plasmon fields amounts only to substitute Y 

_g 

with the running coupling strength, that is ~ (^-^j ^ . Thus, we find that, via renormalization of the 
boundary coupling strength, the collective plasmon modes renormalize the instanton fugacity so as 
to make the gap A between the quantum doublet, | - 7)1,2, and the first excited singlet, |27), scale 

g 

like A(L) = Aq j ^ . Thus, the interaction with the plasmon modes enforces the gap between the 
quantum doublet and the first excited state. 



5.2. Manipulation of the quantum doublet at the FFP 

Let us, firstly, assume that the external phases ju are tuned nearby, but not exactly at, the triple 
degeneracy point, //^/Zj' that is, - iJL\\ln <^ l,|/i2 - fJ-l\/^ ^ 1- The low-energy eff"ective 
Hamiltonian "Kj^ in the basis \R), |G), \B) is now modified to 



ngu 
L 



1 _^ 

9 ^ 4n^ 



.2f^|L -Ye'", _yg-/f 
1+ -Ye^'^ -Ye''^ , (78) 

where we have set 6; = /Ui - fi*,i = 1,2. Of course, for ei = 62 = 0, 'Hrl^] reduces back to the 
Hamiltonian in Eq.(71). In order to rewrite 'Hr\?\ using the states reported in Eq.(74), one has to 
rotate it to Hr\e\ = P^IHrVelP obtaining 



P^'HrVAP 



ngu 



1 f_ 

9 ^ 4n^ 



1 + 



a-i-^ 



V3 



a + i—r= 

V3 



a + i-^ 



27 



V3 



(79) 



-^^,b = -^^z- From Eq.(79), by simply keeping the matrix elements of 



with a = ^-^7^,Ly - -r—r^ 

3 L 2;r' L V3;t 

P^'HfV'^P involving the low energy twofold degenerate ground state, one has 



Eff 



ngu 



1 

9^4^ 



-7U,+ 





a -I- z-7= 

V3 







V3 



(80) 
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From Eq.(80), one sees that, apart from the term proportional to the identity matrix, the matrix 
'Hes can be rewritten as 'Hes = bxcr'' + bycr-', with bx = a, and by = 

To add the z-component to the effective magnetic field one needs to break the degeneracy 
between the | tt)- and the | ||.)-states of So- This may be realized by turning on a nonzero 5|| - see 
Eq.(18)- and amounts to introduce an additional contribution to l-fr given (see Appendix D) by 



B\\ 
4 



z z* 
z* z 
z f 



(81) 



with z = ?e 'f>. When transforming to the states given in Eq.(74), one obtains 



V3 



1 
0-10 




(82) 



From Eq.(82) one readily gets the low energy effective two states Hamiltonian 



'J^Effib] 



ngu 



1 £_ 

9 ^ 4n^ 



-Y\h + b-a 



(83) 



with bx = a, by = and b^ = -^B\\Y. The Hamiltonian 'HEftib] is then the Hamiltonian for 
a spin in an external magnetic field, whose components may be manipulated by acting on the 
external control parameters of the tetrahedral JJN. 

For instance, applying a modulation in time to the phases jl and to the flux /, one may change 
the relative sign between the two states according to the procedure outlined in Ref.[26]. Indeed, 
one may modulate in time 62, so that 62(1) = v sm{a)ot). This results in an effective b field given 
by ^ = (bx,bsm(ojot),bj, with bx,b. constant and b = ^v. The instantaneous eigenvalues of 

'HEsibit)] are then given by +A{t) = + -^Jbl + b\ + lP- %\rL'{(jjQt), while the corresponding adiabatic 
eigenstates are given by 



I - A(0) = e^'-^'^ jcos I - F)i + sin (^J e-'^^l - D^j 

\m) = e^'^^'^ |- sin(^) I - 7). + cos (^j .-'•'^«| - Y)^ , (84) 

with cos(0) = -b^/A{t) and 0(f) = arg(-Z?^ - z^sin(cJoO)- As usual, the phases ^+(t) are chosen 
so as to satisfy the condition 

<A(OI^|A(0) = <-A(OI;^|-A(0) = . (85) 
at at 

Preparing, at f = 0, the system in the state -|A(0)), after twice a period 2T = AnojQ, the relative 
phases between | - y)i and | - 7)2 becomes AO = dt 6(0[1 - cos(0(O)]. Setting b,^ = 
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and V = bx, one finds AO = n and, thus, the relative sign between the two states is exchanged. This 
procedure allows for engineering a NOT port, using the quantum doublet emerging in a tetrahedral 
JJN at the FFP. 

To read the state of the quantum doublet one may look at the pattern of the Josephson currents 
arising circulating in the JJN when it is biased off the triple-degeneracy point. Indeed, within the 
restricted subspace T, the current operators are given by 



^dP^^JM^ e*gu ^ e*gu 



InL 



6L 



n,r = e = ^ ^ + I- 



InL 



3L 



1 1 

1 1 
1 1 

1-1 
-1 1 

1 -1 



(86) 



From Eq.(86), one obtains 



i<-y|/i,^l - y)i = = 2(-y|/i,rl - Y)2 



i<-yiVI-y)2 = 



e gu 
' 6L 



- Y)i = = 2{-Y\hr\ - Y)2 



i{-Y\hr\-Y)2 = i 



e gu 
3L 



(87) 
(88) 
(89) 
(90) 



On the generic state of the doublet given by lor) = cos | |-y)i + e"^ sin 1 1-7)2, the expectation 
value of the current operators are 



/IT I \ ^ / \ ^ 



6L 2 



2nL 



{a\l2,T \a) = — —-{o-^)a + Tr~r^^ 
3L 2 ' 2nL 



(91) 



From Eqs.(B.5) and Eq.(91), it is easy to determine the current pattern identifying each degenerate 
state I -y)i, I -7)2. 



6. Engineering a tetrahedral JJN operating near the FFP 

Spinless TLL leads may be realized also with classical Josephson junctions - i.e. using junc- 
tions for which Ej/Ec > 1- [36], which may be easily and reliably fabricated with well tested 
technologies [32]. In this realization the Luttinger parameter g is given by g ~ yJnEjjEl [36] 
and a tetrahedral JJN operating near the FFP may be fabricated by requiring that g = 2 and then 
setting Lh. ~ 10^; this since, for L > and 1 < g < 3, the phase slip operators destabilize the 
SFP. The requirement g = 2 may be easily satisfied by using junctions for which Ej/E^. ~ 1.3. 
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Choosing A ~ Ej/3 (see Eq.(17)) one gets that L^. ~ y^y ' Lq is of order L^. « 10^ Lq where Lq is 
the reference length of a chain with parameters Ej,Et,. Setting Lq ~ 10 yields then L^, ~ 10^. 

Flux noise induced by the shift {/ ^ / + 5;} affects the stability of the quantum doublet since 
its effects amount to introduce an effective interaction between Sg and a magnetic field B breaking 
the degeneracy between the \R), \G) and |5)-states. As shown in subsection 4, these fluctuations are 
potentially dangerous, as they may induce an effective B-field acting on the spin Sg, which breaks 
the degeneracy between \R), \G) and |5)-states. Since the running coupling strength associated to 
5x scales with L as bj^{L) = LB^ one sees that the scale- Lstop- at which " dangerous" instanton 
trajectories will be suppressed by 5^ may be be defined by the condition Z?x(^stop) ~ j(^stop), from 

3 

which one gets Lstop ~ (^)' Lq. 

To provide a rough estimate of Lstop, one may approximate the actual instanton as a double- 
well instanton, by fitting the parameters of the double-well potential Vdb so that the minimum and 
the maximum points (and the values of Vdb at the corresponding points) coincide with the ones 
obtained from "TYb- A standard computation [23, 37] allows then to estimate the instanton fugacity 
as 

Y -xl exp 

\ gE 

Using the same fabrication parameters as before and fixing ~ 7r/20 (i.e., ~ 0.2£'y(5,)^), one 
gets that Lstop ~ 8 x 10''. As a result, one may infer that noise in the external flux / (described in 
our approach by 5^ and by a finite Lstop) does not affect the quantum doublet provided that 
Lstop > L*. 

7. Concluding Remarks 

In this paper we analyzed the phases accessible to a tetrahedral JJN made by coupling a tetra- 
hedral qubit [26, 27] to three JJ chains acting as TLL leads. 

We showed that, in a pertinent range of the fabrication and control parameters, a robust attrac- 
tive FFP emerges due to the geometry of the tetrahedral JJN . In our approach the central region - 
made by the tetrahedral qubit- is treated as a quantum impurity of this low dimensional network. 
As a result, the emergence of a FFP is a non perturbative phenomenon arising from the strong 
coupling of the impurity with the TLL superconducting leads. 

We argued that the new stable FFP is associated with the emergence of a doubly degenerate 
ground state, which may be regarded as a quantum doublet described by a spin 1 /2 degree of 
freedom, coupled to the plasmon modes of the superconducting TLL leads via the boundary inter- 
action. We showed that this quantum doublet is robust not only against the noise in the external 
control parameters (magnetic flux, gate voltage) but also against the decoherence induced by the 
coupling of the tetrahedral qubit with the superconducting leads. For this purpose, we showed 
that, as the network size increases, the instanton operators, arising from the interaction of the cen- 
tral region with the plasmon modes of the leads, contribute to enforce the energy gap between the 
twofold degenerate ground state and the first excited state. 



-4.36 



3nu 



(92) 
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We also pointed out how one may device protocols allowing to read and manipulate the state 
of the quantum doublet emerging at the FFP; we saw that this may be achieved by connecting the 
tetrahedral JJN to three bulk superconductors at fixed phases {ipj} {j = 1, 2, 3 - see Fig. 2). Indeed, 
we showed that, acting on the {(pj}, induces an "effective magnetic field", which couples to the 
emerging two-level quantum system, providing a tool to prepare the two-level quantum system in 
a given state. 

Finally, it is worth to point out that superconducting devices such as the tetrahedral Josephson 
junction network analyzed in this paper may be used to simulate physical behaviors realizable in 
Kondo systems. Indeed, our RG analysis showed that, forg = 1 and G-(/ = 0) = 0, the tetrahedral 
JJN may be used to simulate a Kondo spin pertinently coupled to band electrons. 
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Appendix A. The central region energy eigenstates 

In this appendix we report the full spectrum of the Hamiltonian given in Eq.(l), for a generic 
value of the applied flux /. In particular, we will single out the twofold degenerate ground state 
whose levels have been used in section 2.2, to define the effective spin- 1/2 operator Sg- The 
eigenstates, together with the corresponding energy eigenvalues, are given by 

m = 2: a fully polarized spin-2 state: 

|2) = imT), (S2 = -2H) (A.l) 

m = 1: four states given by 

11. 1) = , ^ [(-cos / + 1) liTTT) + ITiTT) + ITTiT) + ITTTi)] 

^J2t(t - COS /) 

suif) = -H + ^{-cosf- V3+COS2/) (A.2) 

11. 2) = , ^ [(- cos f-t) liTTT) + ITiTT) + ITTiT) + ITTTi)] 

^J2t{t + cos /) 

su(f) = -// + y (- cos / + V3 + C0SV) (A.3) 



11,3) = -^ [iTiTT) + e'* ITTiT) + e-'"^ ITTTi)] 



'notice that the spin labels correspond to sites 0,1,2,3, respectively; we shall set s - -^T+Vsin^ and f 



-y/3 + cos^ / henceforth; moreover we shall label the energy eigenstates by means of two quantum numbers: the 
former ones refer to the total spin momentum of the states, the latter ones to the z-component of the total spin 
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ei,3(/) = -// + y(cos/- VSsin/) 

11,4) = [lUTT) + e-'"-^ ITTiT) + ITTTi)] 
s,M = -H + ^(cosf+ VSsin/) 



0: six states given by 



10, 1) = 



(liiTT) + e'* liTiT) + e-'* liTTi)) 



•^65(5- VSsin/) 
+ ( V3 sin / - ^) (iTTii) + e'* ITiTi) + e-'* ITiiT)) 



£o,i(/) = Y (cos / - 7l+3sinVj 



10,2) 



yjesis- VSsin/) 

+ (imi) + e-'*iTm) + e'*iTUT)) 

eo,2(/) = eo,i(/) 



( V3 sin/ - 5) (liiTT) + e-'* litit) + e'T litti)) 



10, 3) = 



(liiTT) + e'T liTiT) + e-'T liTTi)) 



-^^65(5 + V3 sin/) 
+ ( V3 sin/ + 5) (iTTii) + e'^ ITiTi) + e-'^ ITiiT)) 

£o,3(/) = y (cos / + ^l+SsinV) 

[( V3 sin/ + (liiTT) + e-'"^ liTiT) + e'"-^ liTTi)) 



10,4) 



-y^6i'(5 + V3 sin/) 
(lTTii) + e-'^ ITiTi) + ITiiT)) 



10,5) = ^ 
V6 



(liiTT) + liTiT) + liTTi)) - (ITTii) + ITiTi) + ITiiT)) 



eo,5(/) = -£7 (cos/ - 1) 

1 10, 6) = ^ [(liiTT) + liTiT) + liTTi)) + (ITTii) + ITiTi) + ITiiT))] 
V6 

eo,6(/) = -£'7(C0S/+ 1) . 



27 



m = -I: four states given by 



1-1, l> = , ^ [(-cos/ + 1) ITiii) + litii) + liiti) + liiiT)] 

^J2t(t-C0Sf) 

£:.u(/) = // + |^(-cos/- V3 + COS2/) (A. 12) 

1-1,2) = , ^ [(-cos/ - 1) ITiii) + liTii) + liiTi) + liiiT)] 
^j2t(t + cosf) 

£_l,2(/) = //+y(-C0S/+ V3+COS2/) (A. 13) 

1-1,3) = ^ [liTii) + e-'* liiTi) + e'T liiiT)] 



V3 

2 



e-i,3(/) = ^ + §^ (cos/ - V3 Sin/) (A.14) 



1-1,4) = -i^ [liTii) + e''^ liiTi) + e-'^ liiiT)] 

e„i,4(/) = ^ + y (cos/ + V3 sin/) . (A.15) 

/n = —2: this is again a fully polarized spin state given by 

|-2) = liiii), (£-2 = 2//) (A.16) 

From the knowledge of these states the effective Hamiltonian for the central region T given in 
section 2. 1 may be easily derived. 

Appendix B. DC-Josephson current pattern at weak couplings 

To induce a DC-Josephson current pattern across the JJN, one has to apply static phase differ- 
ences to the end point of the network. Thus, the currents may be easily computed within imaginary 
time path integral formalism discussed in section 3.1. Indeed, if at the endpoint of branch i a static 
phase (fi is applied, the currents /i, /2, h flowing across the three branches of the JJN may be com- 
puted by taking the logarithmic derivatives of the partition function Z. in Eq.(23) with respect to 
/ii,/Z2 [24], and are given by 



/i = 
h = 
h = 



^ InZ _ 
d\nZ 



dfj.1 
dlnZ 



V6 



- V6 



dlnZ 



dfi2 
d\nZ 
dfi2 



(B.l) 
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To compute Z., one has to sum over the oscillating modes of the fields Xj, by pertinently taking 
into account boundary conditions at both boundaries. At the inner boundary {x = 0), these are 
determined by energy conservation and are given by 



ug dxiO, t) 
2n dx 



--2E, «i sin [aj ■ ^(0, r)] + AE, ^ sin ■ m r) + ^ 



+ 



+ 



AE^Y^j^l 



cos 



+ AE-- 



Yj ^j^G sin 



'2n 

-0--2) 

2n 

-0--2) 



sin[(?y;e(0,T)] + 



(B.2) 



where, in order to account for normal ordering of boundary interaction operators, the boundary 

interaction strengths have been redefined as Ef = (^fY E{, {£ = l,z, 3) [6, 18]. From Eq.(B.2), 
one easily sees that, at the WFP, energy conservation requires Neumann boundary conditions, for 
the plasmon fields at x = (i.e. ^^^^ = ^^^^ = 0). As evidenced in section 4.1, at the outer 
boundary (x = L) xy,X2 obey Dirichlet boundary conditions: ;^'i(L, t) = fXi,X2(L,T) = /i2- As a 
result, the mode expansion of the fields Xj(^^ is 



Xjix,t) =iUj + J- y cos 
fij + (pjix, r) 



n / 1 
1^2'" 



where the oscillator modes «,(«) satisfy the algebra 



[ai(n), ajim)] = 6ij 5„+,„-i,o | " + 2 



(B.3) 



(B.4) 



For our purposes it is convenient to define "spin- 1/2 current operators", that is, operators acting 
on the two-dimensional Hilbert space of So and giving the correct value of the current, when Sg 
is averaged over, as well. It is straightforward to see that these operators are given by 



with 



h = 
h = 
h = 

Zs 



e" I rzd In Zs ^ ird In Zs 



-- V2 ^ 



- V6^ 



^ -V2 



d In Zs 



* 



d/ui 
d lnZs 

dlU2 



V6 



d^2 
dlnZ. 



d^2 



^ ;=i,2 
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(B.5) 



(B.6) 



Relying over the weak coupling assumption, one may sum over the plasmon modes xj within a 
mean-field like approach. In particular, because of the mode expansion in Eq.(B.3), one gets 



<cos[(?, -i^(T)]) =cos[t?, -/l] , <sin[(?, -;^(T)]) = sin[(?, -/l] . (B.7) 

Thus, one readily sees that, resorting to the mean-field approximation amounts to trade x(t) for 
the applied phase differences jl. As a result, one obtains 



-e* din Zs 



2e*Ei ^ aj sin [dj ■ fl] + 4e*E, ^ djS'^ cos [dj ■ fl] 

P J J 

+ 4e*Es cos |y (J - 2) j sin [dj ■ fl] + 

+ 4e*Es sin |y (7 - 2) j sin [dj ■ jl] . 



+ 



(B.8) 



There are two possible ways of interpreting Eq.(B.8): on one hand, one may regard the applied 
phases (and the induced currents) as a probe of the two-level state (which may be set by acting 
upon it with the external fields B\\, 5^). For instance, assuming that the system lies within either 
one of the eigenstates of S^, | H), | H), and computing the average values of the spin operators as 
outlined in Appendix C, from Eq.(B.8) one gets 



< 



d In Zs 
djl 



)s = 2 2 dj{Ei sin \dj ■iJ\+E, cos \dj • /l] } , 



(B.9) 



where the average is computed over the spin coordinates. The corresponding current pattern may 
be derived from Eq.(B.5) and from Eq.(B.9): clearly, it discriminates between the two states | H 
), I U). The same procedure may be applied to probing a generic state of Sg, obviously as long 
as the modification in the state induced by the application of the phases /I is negligible. On the 
other hand, when no other fields are applied to Sg (that is, when 5|| = B^ = 0), the phases (and, of 
course, the currents) themselves may be regarded as defining an effective applied field B, whose 
components are given by 



B 



= -4^3 2 



cos 



B^ = -4^3 sin 



'In 

-jij-2) 
In 

-jU-2) 



cos[ay • yU] 
COS [ay ■ jl\ 



W = +4£, Yj sin[«i ■ fi\ 



(B.IO) 



In this case the phases may be used to drive the state of the two-level system, just as a local 
magnetic field applied to a true spin- 1/2 variable. 
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Appendix C. The imaginary time action for a quantum spin-1/2 variable 

In this appendix we shall review the derivation of the imaginary time path integral formalism 
for a quantum spin-1 /2 variable Sg since it has been used to determine the instanton phases at the 
SFP. 

The starting point is the Euclidean action Z^spin for a spin-1/2 degree of freedom in an external 
magnetic field S, whose dynamics is described by the Hamiltonian 

«Kspin = -S-SG . (C.l) 

A crucial step is the decomposition of the identity I in the basis of the coherent states |0, 0) as 

|<l),0) = e'<^cos(|)|^) + sin(|)|ll) , (C.2) 



The average values of the components of the spin-1/2 variable on the state |0, 0), are then given 
by 



<0, 0|S^|O, 0) = ^ cos(O) sin(0) , <0, 0|S^|O, ©) = ^ sin(O) sin(0) 



1 
2 

In particular, the decomposition of the identity is given by 



<O,0|S^|O,0) = -cos(0) . (C.3) 



1 = J jri 10,0X0, 01 , (C.4) 

with 

dn ... = — J JO J sin(0)J0 ... . (C.5) 

Though coherent states form a complete set, they are not orthogonal to each other; thus, in the 
imaginary time path integral formulation [38], one has to take into account the overlap amplitude 
between two coherent states 

(O. 0, |0„ 0,) = .-'l-'-^l cos cos (I) + sin sin (|) . (C.6) 

Taking into account Eq.(C.6), one gets that the amplitude for Sg to tunnel from the state |Oo, 0o) 
to the state |Oi, 0i) in an (imaginary) time r is given by 



(Oi,0ik-^[^«]|Oo,0o) = J J-U'^'^^Wd-os (e(r)Hj;^/rSX(r)] 

= J £)Q(t) exp[-5£[O,0]] , (C.7) 
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where H[S] is the spin Hamiltonian, while we have defined the polar angles as in Eqs.(C.3), and 
the path integral has to be computed over imaginary time trajectories satisfying the boundary 
conditions 0(0) = 0o, 0(t) = 0i, and 0(0) = Oq, 0(t) = Oi. 

The one- and two-spin imaginary time correlation functions we used in section 3 may be then 
derived by means of a saddle-point approximation. As an example, let us consider the case in 
which the applied field B is uniform and directed along the z-axis, corresponding to Eq.(C.7) for 
the imaginary time amplitude. The saddle-point equations for the functions 0(t), 0(t) are given 
by: 

= ^5|^ = isi„(e(r))[,*(r) + S] 

o0(t) 2 

O = ^^|^ = i0(r)si„(0(r)) . (C.8) 
dO(T) 2 

Eqs.(C.8) imply 

0(t) = constant = 0o = 0i , 0(t) = iBr + Oq . (C.9) 
From Eqs.(C.7,C.9) one finds out that, in the saddle-point approximation, 

<Oi,0i|e-^™|Oo,0o) = 5(00 - 0i)(J(Oy - O,- - iBr) . (C.IO) 
From Eq.(C. 10) and from the identity 



2 /^27r r*7T 

— JOoJOi J0oJ0i sin(0o)sin(0i)<O,,0ik-^™lOo,0o) , (C.ll) 

■^y Jo Jo 



with S^lcr) = cr\cr), one finds out (assuming 5 > 0) that 



{(T\e-'"\o-') = e-"^ . (C.12) 

All the other imaginary time average values listed in Eqs. (32,34a) may be derived following a 
similar approach. For instance, to compute the average value of a component of So, one may use 

<S« (r)) = J](^\e^"Ke-^''\^)e-^^ / E ' ^^-l^) 

O" cr 

and 



1 r^'^ r 

= TT^ d@a d@h sin(0„) sin(0i,) e ^ <cr|0„, 0«)x 

(47r)2 Jo Jo 



<(D„0JS^|«D,,0,)<O,,0,|cr)e^ , (C.14) 

to get the results given in Eqs. (32,34a). When B\\ one may use a similar analysis, provided one 
chooses the z axis directed along the direction of B and rotates the components of Sg accordingly. 
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Appendix D. Derivation of the modulus and phase for instanton trajectories 



In this appendix, we derive the modulus Y and the phase y of the instanton tunnelling ampli- 
tudes. 

First of all we recall that the instanton trajectory may be regarded as the imaginary time evo- 
lution, P = P{t), of the zero mode contribution to Eq.(refeq:mode-exp-strong); it describes a 
tunnelling event between nearest neighboring sites on the lattice of the minima. The "bulk" Eu- 
clidean action for the field x 



4;r Jo Jo 



dx 



u\dT 



dx 



2^ 



(D.l) 



yields 



+ ... = s(°>Dr(T)] + . 



(D.2) 



Wi\hx{T) = x(0,t), M = Lg/(6nu), and Moj^ = ug/{2nL). The ellipses in Eqs.(D.l,D.2) corre- 
sponds to interactions between instantons, mediated by oscillations of the plasmon bulk modes, 
which do not affect the computation of the phase y and will be neglected henceforth. The coupling 
between the ;^-modes and the spin degrees of freedom occurs via the boundary Hamiltonian in 
Eq.(16), which may be presented as 



<K;'^ = S[;e(T)]-SG(T) + 5oDr(T)] 



(D.3) 



with 



7=1 

3 ^2 

!BMr)]=2E,Y,cos -^(j-l) cos[d j ■ j^ir)] 
7=1 L ^ J 

3 r2 



7=1 

3 



(D.4) 



7=1 



Focusing on the "threefold degenerate" point obtained when (//j./i,) = (y.O)^ one easily realizes 
that the instanton trajectories of interest lie along the sides of the triangle whose vertices coincide 
with the R, G, B points defined by nn = «i3 = 0. At the R, G, 5-vertices one gets 



1 



1 



XR = 2n(l,0) , XG = 2n\-, , = 2n\-,— 



V3 



1 1 



2' V3 



(D.5) 
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Any instanton path runs between two of the points in the;^-configuration space listed in Eq.(D.5). 
Moreover, due to 

a2 = 'Ri^]di , a3=5?^(^)«i , (D.6) 



(D.7) 



with 

[-sin(e) cos(0) 

one sees that the path connecting B to R may be obtained by acting with "^(y) on the one con- 
necting G to B; in addition, the path connecting Rto G may be obtained by acting with 'R (^y) on 
the path connecting G toB. As a result, it is enough to compute only one tunnelling amplitude, 
for example the one between G and B. To do so, let us parameterize such an instanton path as 
(Xi('^)^X2(t)) = {n, ^cririj, with cr(0) = and cr(J3) = ^. The corresponding Euclidean action 
is then given by 

S [cr] = dr |^[(o-)2 + (oV^] - 2£i - lE^ ^4-3 sin^ {na)^ . (D.8) 

The term oc cr^ is basically constant, along the instanton path. Thus, the actual value of S [cr] 
may, in principle, be computed by determining the zero-action solution in the "inverted potential", 

O-lnst(T), as 



&'{T)^lE^\^lA-3^m\na,UT)) + l\ , (D.9) 



and, then, evaluating S = 5[crinst]. Since the fugacity F = e ^ is strongly renormalized by the 

g 

interaction with the collective plasmon modes of the bulk, it eventually scales as Y{L) ~ L~3 and, 
thus, one gets that ~ Y(L = Lq). 

The tunnelling amplitudes have also a phase stemming from the topological term due to the 
spin [39]. To derive it, one may consider that the probability amplitude to remain in the same 
state (for example R), is not only given by the instanton/anti-instanton contributions, but also 
by the loops around the three degenerate states. The action of such a loop will induce an extra 
topological term arising from the trajectory of the spin state during the loop. For a generic spin 
state \a(T)) = cos ^ Itl) -i- e"^^'"* sin ^ |JJ.) the topological term can be written as: 

^Top = ^ f dT 0(t) (1 + cos e{T)) . (D.IO) 
^ Jo 

Thus the Euclidean action for a loop will be given by: 

5ioop[;r(T); m; dir)] = Sj^^icp, e] + 35i„,, . (d.i i) 

Writing the spin states along the instanton paths and substituting the instanton coordinates x{t) 
into Eq.(45), one finds the corresponding lowest energy eigenstate. For small one obtains: 

lair)) = cos ^ 1^) + sin ^ 111) (D. 12) 
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where Of is defined in Eq.(5 1) and 7(7) = arg [5jr[;^(T)] + /5j,[;^(t)]J. Since 9 remains constant, the 
evaluation of 5 Top yields 

S Top = in{ 1 + cos Of) ~ in{ 1 - ^) > (D. 1 3) 

that is, the phase contribution to the amplitude of loop tunnelling. Assuming that the three instan- 
ton tunnelling are equivalent, one naturally assign to each tunnelling amplitude a third of the total 
phase. Thus, the single instanton tunnelling amplitude is given by: 

j/ = ye'f-4 , (D.14) 

which is the result used in the paper. 
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